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Abstract 

We give a specific method to solve with quadratic complexity the linear systems 
arising in known algorithms to deal with the sign determination problem. In particular, 
this enable us to improve the complexity bound for sign determination in the univariate 
case. 

1 Introduction 

Let R be a real closed field. A basic problem in computational real algebraic geometry is, 
given a finite set Z C R fc and a finite list V = Pi, . . . , P s of polynomials in R[A 1; . . . , X^], 
to determine the list of sign conditions realized by V on Z. 

A general scheme in most algorithms dealing with the sign determination problem consists on 
the computation of the Tarski-query (also known as Sturm-query) for Z of many products 
of the given polynomials, in order to relate through a linear system this quantities with 
the number of elements in Z satisfying each possible sign condition. For instance, the 
naive approach were the Tarski-query of each of the 3 s polynomials of type Pf 1 . . . Pf s with 
ej = 0, 1, 2 for i = 1, . . . , s is computed, leads to a linear system of size 3 s x 3 s . Nevertheless, 
if m = j^Z at most m sign conditions will be feasible, and then at most m of the coordinates 
of the solution of the considered linear system will be different from 0. 

In [2J, the exponential complexity arising from the number of Tarski-query computations 
and the resolution of the linear system in the approach described above is avoided. This is 
achieved by means of a recursive algorithm where the list V is divided in two sublists, the 
feasible sign conditions for each sublist is computed, and then this information is combined. 
Such combination is obtained by computing at most m 2 Tarski-queries and solving a linear 
system of size at most m 2 x m 2 . 

In [5J, [3] and [U Chapter 10], the methods in [2J are further developed. In [BJ, an algorithm 
is given where the number of points in Z satisfying each feasible sign condition for the list 
Pi, . . . , Pi is computed sequentially for i = 1, . . . , s, following the idea that, at each step, each 
feasible sign condition for Pi, ... , Pj_i may be extended in at most 3 ways. To deal with the 
addition of the polynomial P; to the considered list, at most 2m Tarski-queries are computed 
and a linear system of size at most 3m x 3m is solved. In [3], a more explicit way to choose 
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the polynomials whose Tarski-query is to be computed is given. In [TJ Chapter 10], also the 
feasible sign conditions for Pi on Z are computed at step i, in order to discard beforehand 
some non- feasible sign conditions for Pi, . . . , P, on Z extending feasible sign conditions for 
Pi, . . .,Pi-i on Z. 

Depending on the setting, the Tarski-queries may be computed in different ways, taking 
a different number of operations in the field R, or in a proper domain D containing the 
coefficients of the polynomials in V and polynomials defining the finite set Z. As treated 
with general methods, the linear solving part can be done within 0(m 2 376 ) operations in Q 
(see [I]). In [3j Section 3], where the univariate case is considered, the cost of linear solving 
dominates the overall complexity. In this work, we fix this situation by giving a specific 
method to solve with quadratic complexity the linear systems involved (Theorem Even 
though this method can be used as a subroutine whenever these systems arise, we emphasize 
the result of its use in the univariate case. Following the complexity analisis in [31 Section 
3.3] we obtain: 

Corollary 1 Given P , P x , . . . , P s G R[X], P ^ 0, deg P { < d for i = 0, . . . , s, the feasible 
sign conditions for Pi, . . . , P s on {Pq = 0} (and the number of elements in {Pq = 0} sat- 
isfying each of these sign conditions) can be computed within 0(sd 2 log 3 d) operations in R. 
Moreover, if P has m real roots, this can be done within 0(smd\og(m)\og 2 (d)) operations 
in R. 

The motivation for this work also comes from probabilistic algorithms to determine feasible 
sign conditions in the multivariate case ([5]), which produce a geometric resolution of a set 
of sample points. In this reduction to the univariate case, the degree of the polynomials ob- 
tained equals the Bezout number of some auxiliary polynomial systems, and the complexity 
of the algorithm depends quadratically on this quantity. Using Corollary (U the treatment 
of the univariate case to find the feasible sign conditions for the original multivariate poly- 
nomials can be done without increasing the overall complexity. 

2 Preliminaries and Notation 

We will follow mostly the notation in [H Chapter 10]. In this reference the approach de- 
scribed in Section [T] is followed with the minor difference that the polynomials Pi, . . . ,P S are 
introduced one at each step from back to front; therefore, the notation is adapted to this 
order. 

For i — s, . . . , 1, we call Vi the list Pi, ... , P s , and, at step i, we are given a list E = o~i, . . . , o~ r 
of elements in {0, 1, — l} Vi with oi <i cx • ■ • <i cx a r (0 -< 1 -< — 1) containing, may be properly, 
all the feasible sign condition for Vi on Z and we are to compute the exact list of feasible 
sign condition for Vi on Z and the number of elements in Z satisfying each of these sign 
conditions. Note that the inequality r < 3m holds at every step. 
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For P G R[X], we note c(P = 0, Z), c(P > 0, Z) and c{P < 0, Z) the number of elements in 
Z satisfying the condition P = 0, P > and P < respectively. Recall that the Tarski-query 
of a polynomial P for Z is the number 

TaQ(P, Z) = c(P >0,Z)- c(P <0,Z). 

For a G {0, 1, we note c(a, Z) the number of elements x in Z satisfying sign(Pj(x)) = 

Oj for j = i,...,s and, for a list £' = o' x ,...,a[ of elements in {0,1,-1}^, we note 
c(H',Z) the vector whose components are c(a[, Z), . . . , c{a'^ Z). We also note a the ele- 
ment of {0, 1, — l} Vi+1 obtained from a by deleting the coordinate corresponding to Pi and 
E' the list a[, . . . , a[. 

We divide the given list E in 12 ordered sublists as follows: for ^ B = {b±, . . . , bi) C 
{0, 1, —1} with bi -< • • • -< bi and for b G B, the list E[| b is composed by those o G E such 
that cr(Pj) = 6 and the set 

{b' G {0, 1, —1} | 3a' G E such that a' extends a and cx'(Pj) = b'} 

equals B. For simplicity we note E , Ei and E_i for E[j, Ej and El x respectively. In addition, 
since E^ ^ is the same list for every b G B, we note E^...^ for any such list. 

We also divide the list E in 3 ordered sublists as follows: E(i) is the merge of lists E , Si, £_i, 
Eq 1; Eq_ 1; E^lj and Eq^, E( 2 ) is the merge of lists Eq 1; T.q^, Sj_ 1 and E^.]^ and E( 3 ) is 
the same list than Eqj 1 ^. 

Consider also the list Ada(E) of elements in {0, 1,2}''* (which represents a list of multide- 
grees) defined recursively by: 

if i = s and r = 1, 

0, 1 if i = s and r = 2, 

0, 1, 2 if i = s and r = 3, 

x Ada(E (1) ), 1 x Ada(E (2 )), 2 x Ada(E (3) ) if i < s. 

For a list A = ax, . . . , of elements in {0, 1, 2} Vi and a list E' = a[, . . . , o[ 2 of elements in 
{0, 1, -l} n , we note Mat(A, £') the Z hxh matrix defined by 

Mat(AE') ili2 =<7^\ 
for ji = 1, . . . , li and j 2 = 1, . . . , l 2 , with the understanding that 0° = 1. 



Remark 2 Following this notation, we have that: 

Mat(Ada(E), E) c(E, Z) = TaQ(P, Ada(E) , Z) (1) 

and the rxr matrix Mat (Ada(E), E) is invertible (see |3J Sections 2 and 3} or |H Proposition 
10.65]). 
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For a matrix M with rows indexed by a list A of multidegrees and columns indexed by a list 
£' of sign conditions, and for any sublists A' of A and S" of E', we will note, if convenient, 
Ma'i My," and Ma>,h» the submatrices obtained from M by taking only the rows in A', only 
the columns in E", and only the rows in A' and the columns in E" respectively. We will use 
a similar notation for vectors whose coordinates are indexed by a list of multidegrees or a 
list of sign conditions. 



3 The specific method for linear solving 

Note that a different order in E would lead to a permutation of columns in the matrix 
Mat(Ada(E), E) and the elements of the vector c(E, Z). To explain our method in a simpler 
way, we will suppose that we change the order in E in such a way that we find first the 
elements in E(i), then those in E( 2 ) and finally those in E( 3 ). Nevertheless, this change of 
order is not actually necessary in the execution of the linear solving algorithm. 

If % — s, we have that either r = 1; r = 2 and E = 0, 1; r = 2 and E = 0, —1; r = 2 and 
E = 1, — 1 or r = 3. Depending on which of these conditions holds, Mat(Ada(E), E) is one 
of the following matrices: 




If i < s, consider the matrices Mi = Mat(Ada(£(i)), £(i)), M 2 = Mat(Ada(E( 2 )), E( 2 )) and 
M 3 = Mat(Ada(£ (3) ),£ (3) ). Then, Mat(Ada(E), E) is the matrix: 



Mi 


M[ 


M'l 


X 


M 2 


-M> 2 


Y 


Z 


M 3 



\ 



where: 



• M[ is the matrix formed with the columns of M x corresponding to sign conditions in 
Eoi, Eo-i, Ei_i and Eoi-i, 

• M" is the matrix formed with the columns of Mi corresponding to sign conditions in 
Eoi-i, 

• M 2 is the matrix obtained from M 2 multiplying by —1 the columns corresponding to 
sign conditions in £ _i, 
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• M' 2 is the matrix formed with the columns of M2 corresponding to sign conditions in 

E01-1, 

• X = Mat(l x Ada(S( 2 )), Em), Y = Mat(2 x Ada(E (3) ), Em) and Z = Mat(2 x 
Ada(E (3) ), E( 2 )). 

Remark 3 • The only non-zero columns in matrices X and Y are those corresponding 
to sign conditions in Si, E_i and E^. 

• The following relations are satisfied: 



• Since E(3) is included in £(2), Ada(E( 3 )) is included in Ada(E(2)) and then we have: 

^lxAda(S (3) ),Si = ^Sl! -**lxAda(E (3) ),£-i = — ^£-i> ^lxAda(S (3) ),E~i 1 = — ^-i' 

In order to describe our method to solve linear system (TO , we define the following matrices: 
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where 1%, I2 and ^3 denote the identity matrices which size is the length of Em, E( 2 ) and E( 3 ) 
respectively and, if we index the columns of Ji, J 2 and I 3 with Em, E( 2 ) and E( 3 ): 
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• 1% is the matrix obtained from I2 multiplying by —1 the columns corresponding to sign 
conditions in £0-1 and by | the columns corresponding to sign conditions in £i_i, 

• Z is the matrix obtained from Z multiplying by the columns corresponding to sign 
conditions in E}_ 1; 

• I' 2 is the matrix formed with the columns of I2 corresponding to sign conditions in 
S01-1, 

• I[ is the matrix formed with the columns of I\ corresponding to sign conditions in 
Sqi, E0-1, and £oi-i> 

• is the matrix formed with the columns of l\ corresponding to sign conditions in 

Soi-l- 

Proposition 4 The matrix Mat(Ada(£), £) _1 equals the product Ng . . . N\. 

Proof: First, note that 

Mf 1 M[ = I[, M^ l M'( = M 2 ~ 1 M' 2 = I 2 . 

Because of the first item of Remark [3] we can conclude that 

Xl'{ = 0, YI'{ = 0, 

and using the first and second items of Remark [3], we can conclude that 

-XI[ + M 2 = M 2 , -YI[ + Z = Z, ZI' 2 = M 3 , 

where M 2 is the matrix obtained from M 2 multiplying by —1 the columns correspond- 
ing to sign conditions in £0-1 an d by 2 the columns corresponding to sign conditions in 
Ei_i. With all these relations, the proof can be done by simple computation of the product 
Ng... iViMat(Ada(£), £). 

□ 

The decomposition of Mat(Ada(S), S) _1 as a product of several matrices in Proposition H] 
leads us to the following recursive algorithm: 



Algorithm: Auxlinsolve 
Input: S,TaQ(pf da(E) ,Z). 
Output: c(S, Z). 
Procedure: 
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1. If z = s, return Mat(Ada(E), £)- 1 TaQ(7^ da(E) , Z). 

2. If i < s: 

0. Initialize c = TaQ(pf da(S) , Z). 

1. c S{1) = Auxlinsolve(E(i),c S(1) ). 

2- c S{2) = c E(2) -Mat(l x Ada(E (2) ), £i)c Sl -Mat(l x Ada(E (2) ), E^c^ -Mat(l x 
Ada(S (2) ),S 1 l 1 1 )c sr _i i ; 

cs (3) = c S(3) -Mat(2 x Ada(E (3) ), E x )c Sl - Mat (2 x Ada(E (3) ), E_i)c s _ 1 -Mat (2 x 
Ada(E (3) ),E^! 1 )c E -_i i . 

3. c S{2) = Auxlinsolve(E( 2 ), c S(2) ). 

4. c v -i = — c r -i ; 

^0-1 ^0-1 

5. c S(3) = c S(3) - Mat(2 x Ada(E (3) ), S^c^ - Mat (2 x Ada(E (3) ), E^c^ - 
M 3 c s i . 

6. c S(3) = Auxlinsolve(E (3) ,c E(3) ). 

^' CS (3) ~~ 2 CE (3)- 

8. Cyl = Cyl + C E ,_, . 

9. CyO = CyO — Cyl \ 

^01 ^01 ^01 ' 

Cyo = c^o — c v -i; 

^1-1 ^1-1' 

CyO = CyO — Cyl — Cv„, . 
^01-1 ^01-1 -^01-1 ^(3) 

10. Return c. 



Theorem 5 Algorithm Auxilinsolve solves linear system (QJ) within 2r 2 operations in Q. 

Proof: The correctness of the algorithm follows since, for j — 0, ... ,9, after Step 2.j we have 
computed 

c = Nj ... AT 1 Mat(Ada(E), E)c(E, Z). 

To bound the number of operations needed, we proceed by induction on i. If i = s, the 
result follows, since the inverse of the 5 possible matrices Mat(Ada(E), E) is pre-computed 
and the product by TaQ(Pi Vda( ' E ' ) , Z) takes r(2r — 1) operations in Q. 

Suppose now i < s. For 7^ £> = {61, . . . , b{\ C {0, 1, —1} note by r^...^ the size of the list 
E^ 6 for any b G B. Using the inductive hypothesis, the number of operations in each step 
is bounded in the following way: 
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2.1. 2(r + ri + r_i + r i + r _i + n_i + r i_i) 2 . 

2.2. 2(r i + r _i + r x _i + 2r Q1 ^ l )(r 1 + r_i + n_i). 

2.3. 2(r i + r _i + n_i + r i_i) 2 . 

2.4. ro-i+ri_i. 

2.5. 2r i_i(r i + r _i + r i_i). 

2.6. 2r 2 1 _ 1 . 

2.7. roi-i. 

2.8. roi-i- 

2.9. r i + r _i + n_ x + 2r i-i. 

Since the sum of all this numbers is always lower than or equal to 2r 2 = 2(ro + r\ + r_i + 
2r i + 2r -i + 2ri_i + 3r i_i) 2 , the result follows. 

□ 

Remark 6 The third item of Remark^ implies that Step 2.2 could be replaced in the fol- 
lowing way: 

2.2.' v = Mat(l x Ada(E( 2) ),Si)c Sl ; 
v' = Mat(l x Ada(S (2) ), S_i)c s _ 1 ; 
v" = Mat(l x Ada(E( 2 )),Sr^ 1 )c s -i i ; 
c S(2) = c S(2) -v-v'- v"; 

Cs (3) = C£ (3) — flxAda(E (3) ) + t, lxAda(E (3) ) + V lxAda(E (3) )- 

which takes a smaller number of operations than Step 2.2. 
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